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arstpaot  » shortcomings  of  the  n/d  method  and  of  the  Chew-Low-Salzman 
method  for  the  numerical  solution  of  singular  nonlinear  integral  equations 
of  the  Low  type,  utilizing  the  preliminary  regularization  of  those  equations, 
are  analyzed.  The  method  is  presented  for  the  direct  numerical  solution 
(without  preliminary  regularization)  of  one  particular  class  of  Low  type 
integral  equatlons--the  Chew-Low  equation,  which  describes  the  p-wave 
resulting  from  the  pion-nucleon  interaction  under  the  assumption  that  the 
nucleon  is  fixed.  By  utilizing  certain  substitutions  of  variables,  the 
Chew-Low  equation  Is  transformed  to  a  form  for  which  solution*.  are  derived 
by  the  method  of  successive  approximations.  Measures  for  eliminating  the 
instability  from  the  Iterative  process  are  indicated.  The  authors  explain 
how  the  proposed  method  can  be  used  to  obtain  adiabatic  and  resonance 
solutions  of  the  Chew-Low  equation.  It  is  shown  that  the  adiabatic  solution 
ceases  to  be  analytic  when  the  coupling  constant  f2  =  O.07.  Some  numerical 
results  obtained  by  authors  at  the  Computing  Center  of  the  Joint  Institute 
of  Nuclear  Research  at  Dubna  are  presented. 

Orig.  art.  has:  2  figures  and  11  formulas.:  English  translation:  9  pages. 
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THE  NUMERICAL  SOLUTION  OF  ONE  CLASS  OF  NONLINEAR  INTEGRAL 
EQUATIONS  OF  DISPERSED  TYPE 


I.  Nedyatkov  and  G.  Penchev 


This  work  presents  one  method  of  numerical  solution  of  nonregularized 
Low  integral  equations.  Low  equations  are  systems  of  ^nonlinear 
integral  equations  of  dispersed  type  (which  include  Cauchy  integrals). 

For  various  values  of  N,  particular  cases  of  Low  equations  obtain. 
Thus. for  example,  where  N  is  equal  to  1  or  2,  we  have  Castile jo,  Dalitz 
and  Dyson  equations,  where  N»2  --  Wanders  equations,  where  N*3,  Chew  and 
Low  equations  [1],  Chew  and  Mandelstam  equations  [3],  Shirkov  equations 
[5],  etc. 

In  modern  investigations  of  the  numerical  solution  of  Low  equations, 
two  methods  are  used:  the  ND-method  and  the  similar  Chew -Low-Sal zaan 
method  (2].  It  is  characteristic  for  both  methods  that  first  the  Low 
equations  are  regularized,  then  solved  by  iteration.  An  algebraic  method 
essentially  different  from  these  two  was  recently  suggested,  but  we 
will  not  discuss  it  [6]. 

It  is  known  from  the  theory  of  singular  integral  equations  [8] 
that  the  solution  to  regularized  equations  may  not  correspond  to  the 
solution  to  the  Initial  equations  and  due  to  this, in  each  individual 
case .detailed  investigation  is  required.  This  is  also  true  of  the  two 
methods  mentioned  above.  Due  to  the  fact  that  the  denominators  Du(z)  and 

ha(z)  may  have  zeros  in  the  complex  equation  sometimes  solutions  are 
produced  which  do  not  satisfy  the  original  equations.  Since  in  practical 
calculations  the  zeros  in  the  denominator  are  not  sought  out,  the  usage 
of  these  two  methods  involves  a  certain  risk. 

The  purpose  of  the  present  work  is  to  present  a  numerical  method 
for  solution  of  Low  equations  which  does  not  have  the  defects  of  the 
ND-method  or  the  Chew-Low-Salzman  method.  Just  as  in  the  case  of  linear 
singular  integral  equations  [9,  10,  11],  the  procedure  involves  direct 
solution  of  the  initial  equations  without  preliminary  regularization. 

Due  to  the  presence  of  nonlinearities  in  the  integral  equations,  the 
difficulties  of  numerical  solution  are  much  greater  than  those  in  (9,  10, 
11]. 


Since  the  difference  in  the  individual  form  of  the  equations  is 
inessential  from  the  point  of  view  of  the  method,  throughout  the  following 
we  will  concern  ourselves  with  a  single  special  type  of  Low  equation  --  the 
Chew-Low  equation  [1],  using  which  one  can  describe  the  P-wave  of 
pion-nucleon  scattering  on  the  assumption  that  the  nucleon  is  fixed. 

In  section  1  below  we  give  formulas  for  replacement  of  variables 
and  outline  the  primary  operational  formulas ;  in  sevens  2  and  3  we 
present  the  necessary  explanation  concerning  the  application  of  the  method 
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of  sequential  approximations  to  produce  adiabatic  and  resonant  solutions, 
and  in  section  4  a  few  numerical  results  which  we  produced  using  an 
electronic  conputer  in  the  computer  center  of  the  Joint  Institute  of 
Nuclear  Research  in  Dubna  are  presented. 


1.  Conversion  of  the  Equation 

The  Chew- Low  equation  can  be  written  in  the  form 


1  fv  >  v  ['  W  i* ,  v  A i  1  , 


(l) 


where  h  are  unknown  functions,  analytic  in  plane  z  with  cross-section 

az 

(-•»,  -1)  and  (-fl,  +«)  with  the  exception  of  point  z*0,  where  they  have 
simple  poles  with  residuum  A  ;  z-*w+i e  where  w>l,  z+m-ie  where  w<-l, 
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p.fw2-!)1/2  is  the  impulse  of  the  pion  in  the  system  at  the  inertial  center, 
u  is  the  total  energy  of  the  pion  in  the  same  system,  f2  is  the  coupling 

2  2 

factor,  v(p)«l/(l+p  a  )  is  the  formation  function  (a*0.27)  with  the  given 
h,  c  and  rest  mass  of  the  pion  in  the  position  equal  to  unity,  and  index  a 
has  the  value 


a*  I  — * 
a«=2— ♦ 
«^-3  — * 


The  first  index  in  the  parentheses  is  the  total  isotropic  spin  of 
the  pion-nucleon  system,  the  second  index  --  is  the  total  moment  of  the 
system. 

From  (1)  and  (2),  certain  properties  of  function  h  (u>)  follow 

EL 

immediately: 
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This  can  be  written  as  follows: 
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in  which 
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(9) 


If  we  represent  this  portion  of  equation  (8a)  as  Q  (t) ,  obviously 


m±)rL\t)-4Ql(() 

2 - 


(8b) 


Finally,  multiplying  (7)  by  tg  6a (t)  and  using  the  representation  of 
(9) ,  we  produce 


(10) 
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Functions  h^fu)  are  connected  with  the  phases  of  scattering  dft(u)  of 
the  P«wave  by  relationship 


cxp|/^-OlHiu^H 
- >•»(*)  o»(,„)  - 


r»»  -  J . 


(4) 


To  find  functions  R  h  (u)  and  J_h  (u»)  (u»>l)  it  is  necessary  to 

“  A  m  g 

solve  equations  (3a)  and  (3b)  simultaneously,  which  produces  a  system  of 
three  nonlinear  singular  integral  equations  and  three  algebraic  equations . 
We  will  use  the  relationship  (4): 


Sill  M<'>)  COS  4.(#i)  in  .  I  dot  SflJ*  6m((>t) 

.  +  *  ^ 
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Now  let  us  introduce  the  boundaries  of  integration  by  replacement  of 
variables 


1 

l 


(6) 


For  simplicity  after  the  replacement  of  variables, we  introduced  new 
symbols  for  the  functions  h&,  6^,  v  and  p. 

Substituting  (6)  in  (5) ,  we  produce  the  system 


-  3  - 
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In  working  for  an  adiabatic  solution,  the  sign  before  the  square  root  in 
(8b)  oust  be  negative,  and  in  the  third  phase  (a*3)  of  the  resonant 
solution  it  was  proven  that  beginning  with  the  minus  on  even  numbered 
paths  the  same  sign  changes  [4].  In  the  case  which  we  have  analyzed,  the 
sign  of  two  paths  changed. 

In  [5]  it  is  shown  that  in  another  type  of  equation  (the  equation 
describing  pion-pion  scattering)  three  types  of  solutions  are  possible  for 
various  asymptotic  behaviors  with  increasing  energy  (1/ln  u,  l/u>  and 

2 

l/u >  ).  We  will  seek  the  solution  which  approaches  zero  at  infinity. 


2.  Adiabatic  Solution 

We  note  the  fniinwit*®  fact.  It  is  obvious  that  sys teas  (8)  and  (10) 
are  equivalent,  i.e.  that  if  calculations  are  performed  with  sufficiently 
hitfi  accuracy,  the  solutions  produced  will  be  the  same  regardless  of  the 
system  used.  However,  due  to  limitations  on  the  accuracy,  the  systems 
may  be  found  inadequately  effective.  We  believe,  for  example,  that  the 
accuracy  with  which  calculations  are  performed  in  this.work  is  insuffi¬ 
cient  for  an  adiabatic  solution  to  the  problem  where  ^>0.05,  since  near 
the  maximum  of  function  X^(t)  the  expression  under  the  radical  in  (8b) 

becomes  negative.  Due  to  this,  we  will  seek  the  adiabatic  solution  by 
solving  system  (10)  and  using  the  first  relationship  in  (9).  To  use  the 
method  of  sequential  approximations  we  must  know  the  coupling  factor 

f2  and  use  the  first  approximation  of  asplitude  fig  (*>1,2,5).  Since  the 

speed  of  convergence  of  the  iterational  process  depends  on  the  suitability 
of  the  first  approximation,  usually  a  program  is  written  which,  even 
though  the  convergence  may  not  be  too  rapid,  has  advantages  as  concern* 
systematic  ordering.  This  sort  of  program  is  required  in  correspondence 
with  ,the  well  known  concept  of  the  adiabatic  interaction  in  which  at 

f2-*)  all  phase  differences  are  equal  to  zero.  Then  the  first  approximation 

with  f2  can  be  produced  using  (8b),  assuming  The  subsequent 

approximations  with  the  f2  selected  are  produced  by  solving  system  (10). 

After  this,  varying  f2  with,  sufficiently  small  increment  Af2  with  the  value 

of  f2  as  the  zero  approximation,  we  use  the  solutions  produced  corresponding 

to  tiia  preceding  f2 . 


I.  Meson  ant  Solution 

Our  calculations  have  shown  that  the  adiabatic  solution  cassot  be 
coordinated  with  experimental  results  even  qualitatively  as  calculated  in 
this  work;  no  therefore  sot  ourselves  the  task  of  seeking  the  resonant 
solution.  However,  since  there  can  he  m  unlimited  mmd>er  of  »uch 
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solutions,  we  sought  the  primary  resonant  solution,  corresponding 
closely  to  the  experimental  curve.  This  means  that  we  must  find  a 
solution  in  which  the  first  two  phases  are  sr  11  in  comparison  with  the 
third  phase,  in  which  the  resonance  is  located.  The  program  which  we  used 
in  producing  the  adiabatic  solution  is  not  suitable  for  producing  the 
resonant  solution.  With  the  given  values  of  coupling  factor  as  the  first 
approximation  of  solitudes  6i  and  62  we  can  use  6j(t)=0  and  d^tUO. 

The  first  approximation  of  the  third  phase  is  produced  from 


Ctgda(a»)R: 


3ftf*(w0— of) 

vK(pf 


(ID 


where  wq-2.17,  f^g«<0.Q87,  M*6.67  (mass  of  the  nucleon)  and 

Hits  last  formula,  as  is  shown  in  the  work  by  Layson  [7],  gives 
good  correspondence  with  the  results  produced  for  experimentally. 

As  was  mentioned  in  section  2,  systems  (8)  and  (10)  may  be  found 
inadequately  effective.  This  appears  particularly  strong  in  an  attempt 

to  find  the  resonant  solution.  The  numerical  investigations  which  we 
performed  an  our  electronic  cosputer  showed  that  it  is  impossible  by 
iterat  »  to  produce  the  resonant  solution  using  equations  (8a)  and  (8b). 

We  have  already  noted  that  the  first  two  phases  require  the  minus 
sign  before  th.  spn  root,  nod  th.t  iy  duo  to  th.  ro.on.nce.  wquires 

two  paths  with  variable  sign.  It  was  found  that  the  points  where  the 
sign  changes  (where  t*0.4  and  t»0.57)  discontinuities  appear  which 
disrupt  the  smooth  course  of  the  curve.  Equation  (10)  is  essentially 
inapplicable  to  tine  resonant  solution,  since  it  does  not  contain  tg  6  , 

and  at  t«0.5(w»2)d j**/2.  Due  to  this,  we  decided  to  use  a  combined  method 

which  consisted  of  the  following:  at  the  beginning  of  the  interval,  up 
to  t«0,4S,  the  calculation  was  conducted  according  to  (10)  and  the  first 
formula  in  9).  then  in  the  interval  from  t«0.4S  to  t«0*SS  --  according 
to  (la),  (ft>)  and  (9)*  in  which  in  (lb)  the  minus  sign  was  used  for 
5 j  and  ty  the  pirns  sign  for  d^;  in  the  next  interval  calculations  were 

performed  once  awe  according  to  (10)  nod  the  first  formula  of  (9). 
this  was  dome  in  order  to  eliminate  the  deficiencies  in  systems  (8)  and 
liu; .  la  the  calculation  process  it  was  fotnd  that  it  was  necessary  to 
artificially  snoothe  function  dj  in  those  places  where  the  equations 

were  exchanged,  since  with  increasing  number  of  iterations  at  these  points 
(t«0.4S  and  t*0.SS)  discontinuities  gradually  appear. 


4.  Hum r lea)  Hmselts 

Calculations  men  conducted  using  an  electronic  computer  and  the 
methods  described  in  sections  1,  2  and  S.  The  interval  from  0  to  1  for  t 
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was  divided  into  100  sub intervals.  Also,  to  improve  the  accuracy,  the 
integrals  were  calculated  with  automatic  selection  of  the  step;  at  points 
other  than  those  selected,  interpolation  was  performed.  The  integrals 
were  calculated  with  a  relative  accuracy  of  S' 10  ,  and  iteration  continued 

until  the  modulus  of  the  difference  between  two  successive  iterations  was 

less  than  5 '10  .  Approximately  five  minutes  of  machine  time  was 

required  for  each  complete  cycle. 


GRAPHIC  NOT 


Figure  1. 


Figure  2. 


In  the  calculation  of  the  adiabatic  solution,  the  coupling  factor 
varied  by  Af*«0.01.  The  iterations  converged  where  f*<0.07.  In  this  case, 

when  f2  increased,  the  nuatoer  of  iterations  also  increased,  increasing 

froa  seven  iterations  at  f^-O.Ol  to  21  iterations  at  f^-O.O?.  be  aust 

alto  note  that  where  f**Q.07  the  following  phono— non  was  observed:  tp  to 
to  the  thirteenth  iteration  the  process  converged  noraally,  while  between 
the  fourteenth  aad  sixteenth  iterations  two  peaks  appeared  in  the  function 
(at  t-0.14  aad  t»0.22)  after  which  they  converged  to  a  single  peak  at 
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t=0.18  and  the  result  after  the  twenty-second  iteration  approached  the 
solution  of  the  problem.  When  the  iterational  process  was  continued 

beyond  f2.*  .07,  the  same  phenomenon  was  repeated  and  as  was  to  be 

expected  at  f^=0.08  the  sequential  approximations  no  longer  converged. 

All  of  the  functions  produced  and  subjected  to  checking  approximately 
satisfied  the  initial  equation.  It  was  found  that  in  the  adiabatic  case 
the  solutions  which  we  determined  satisfied  crossing  symmetry  (3)  with  an 
accuracy  of  0.011,  while  the  resonance  solution  showed  an  accuracy  of 
about  O.H.  We  present  below  some  of  our  solutions  for  various  values 
of  coupling  factor.  Figure  1  shows  S  (w) (a»l,2,3)  in  degrees  for  the 

adiabatic  case,  figure  2  —  6 ^ (oj)  and  Sjta)  for  the  resonant  solution; 

<S3(w)  for  the  same  solution  is  shown  in  table  1. 

The  authors  express  their  gratitude  to  I.  Todorov  and  R.  Genshev 
for  useful  discussion. 


Table  1 
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DISTRIBUTION  LIST 


Organization  Nr .  Cys 


AIR  FORCE 
Hq  USAF 

ACIC  (ACDEL-7)  2 

OAR  (RRYA)  1 

AFCRL  (CRXLR)  1 

ARL  (ARB)  WP.AFB  1 
Hq  AFSC  (SOFT)  1 

AEDC  (AEY)  1 

AFWL  (WLF)  1 

ASD  (ASFS  -2)  8 

SAMSO  (SMFA)  1 

FTD 

TDBAE  3 

ATD  (2) 

PHS  (1) 

TDBAS-2  2 

TDBID-2  2 

TDBR  1 

TDDTE  1 

TDGS  1 


TDP  (PHE/PRQ  1 


Organization  Nr.  Cya. 


OTHER  POD  AGENCIES 

DX  ?-0 
Hq  USAREUR  1 
Army  Map  Service  1 
Msl  Intel  Dir  RDST  1 
U.S.  Army  (FSTC)  * 
Army  Security  Agency  1 
Harry  Diamond  Lab  1 
NOTS  China  Lake  1 
PAC  Msl  Range  1 


OTHER  GOVERNMENT  AGENCIES 

AEC  (Tenn)  2 

AEC  (Wash)  2 

NAFEC  1 

NASA  (ATSS-T)  1 


DISTRIBUTION  TO  BE  MADE  BY 


DIA  (DIACO-3) 

DIAAP-1C1  (1) 

DIAAP-1H2  (1) 

DIAST-A  U) 

DIAST-5C  (1) 

U.S.  Navy  (ONR)  (1) 

0P922F2  (1) 

CIA  (SD)  (5) 

NS  A  (CREF/CDB)  (6) 
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